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ABSTRACT 

The use of large, X-ray selected galaxy cluster catalogues for cosmological analyses 
requires a thorough understanding of the X-ray mass estimates.Weak gravitational 
leasing is an ideal method to shed light on such issues, due to its insensitivity to the 
cluster dynamical state. We perform a weak leasing calibration of 166 galaxy clusters 
from the REFLEX and BCS cluster catalogue and compare our results to the X-ray 
masses based on scaled luminosities from that catalogue. To interpret the weak leasing 
signal in terms of cluster masses, we compare the leasing signal to simple theoretical 
Navarro-Frenk-White models and to simulated cluster leasing profiles, including com¬ 
plications such as cluster substructure, projected large-scale structure, and Eddington 
bias. We find evidence of underestimation in the X-ray masses, as expected, with 
(Mx/Mwl) = O.TS^o^os best-fit model. The biases in cosmological parameters 

in a typical cluster abundance measurement that ignores this mass bias will typically 
exceed the statistical errors. 

Key words: galaxies: clusters: general - large-scale structure of Universe — methods: 
numerical — gravitational leasing: weak — X-rays: galaxies: clusters 


1 INTRODUCTION 

Galaxy clusters are the most massive objects in the universe, 
reaching masses of or more. Their abundances are 

a sensitive probe of the growth of structure because they 
measure the exponential tail of the halo mass function (e.g. 
Bahcall 1988; Vikhlinin et al. 2009a). They also provide an 
important cosmological probe, as they are complementary 
to the traditional geometric probes (Albrecht et al. 2006) 
and can be used to test modified theories of gravity (Jain 
& Zhang 2008). However, for robust interpretation of clus¬ 
ter abundance measurements we must know the masses of 
the clusters in the sample; uncertainties in inferring clus¬ 
ter masses from observations are a major limitation to the 
usefulness of cluster abundances in precision cosmological 
measurements (Vikhlinin et al. 2009b; Mandelbaum et al. 
2010; Allen et al. 2011; Hoekstra et al. 2013; Weinberg et al. 
2013; von der Linden et al. 2014a). 

Cluster mass measurements or mass proxies can be de¬ 
fined from measurements in multiple wavebands, deriving 
from different components of the cluster structure. From 
simulation and observation, we know that the dominant 
component of galaxy clusters is dark matter, usually dis¬ 
tributed in a somewhat clumpy, roughly triaxial structure 
with a density profile well-described by a Navarro-Frenk- 


White profile (NFW; Navarro et al. 1997) or the related 
Einasto profile (Navarro et al. 2010). The dominant bary- 
onic component of clusters is a reservoir of hot, ionized gas 
that collects in the potential well of the cluster. In the X- 
ray, this gas is clearly visible as a bright, extended extra- 
galactic source. In the radio, the gas also appears as devi¬ 
ations in the cosmic microwave background, known as the 
Sunyaev-Zel’dovich effect (Sunyaev & Zeldovich 1970): the 
hot gas Compton upscatters CMB photons, leaving a decre¬ 
ment or increment in the CMB brightness depending on 
wavelength. Finally, in the optical, galaxy clusters are char¬ 
acterized by an overdensity of galaxies in a small region of 
physical and redshift space. Mass estimates can be made 
using any of these wavebands as proxies: scalings to optical 
richness (number of galaxies in the overdensity), strength of 
the Sunyaev-Zel’dovich effect, and X-ray quantities such as 
temperature and luminosity have all been used, in addition 
to the weak lensing measurements discussed below. 

Numerical studies indicate that X-ray masses which as¬ 
sume hydrostatic equilibrium are biased low at the 10-30 
per cent level depending on assumptions about feedback and 
preheating (e.g., Rasia et al. 2006; Lau et al. 2009; Battaglia 
et al. 2012; Rasia et al. 2012; Kay et al. 2012; Nelson et al. 
2014, and references therein). This is generally called the 
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“hydrostatic mass bias.” In addition, errors may be intro¬ 
duced through differing X-ray normalizations and measure¬ 
ments. The total offset between X-ray mass measurements 
or mass scalings and weak lensing masses is usually mea¬ 
sured as (Mx/Mwl) — 1 or as the slope of an Mwl vs Mx 
relation minus one; we will call these two cases the “average 
bias” and the “fit bias.” 

Several groups have found these mass biases to be dif¬ 
ferent for so-called “relaxed” and “unrelaxed”, or “undis¬ 
turbed” and “disturbed”, clusters. The division is based on 
the amount of substructure present in the X-ray gas maps: 
clusters with little substructure are relaxed or undisturbed, 
those with a high degree of substructure are unrelaxed or 
disturbed. Unrelaxed clusters are thought to indicate recent 
mergers. The cluster masses in X-ray are determined as hy¬ 
drostatic masses (a combination of the X-ray temperature 
and gas or electron densities out to a given aperture, as¬ 
suming the gas is fully in hydrostatic equilibrium and often 
that it is spherically symmetric; see e.g. Israel et al. 2014 for 
details) or by a calibrated relationship between total mass 
and a quantity such as luminosity or temperature, which we 
will call “mass scalings” in this work to distinguish them 
from mass measurements made using multiple observational 
quantities. The aperture is defined by the overdensity A, 
the ratio by which the average density inside the aperture 
is greater than the critical density pc or (more rarely) the 
mean matter density pm- For clarity these are often called 
Ac and Am, respectively. Common choices are A = 500 for 
X-ray analyses (which typically cannot detect gas very far 
out from the center of the cluster) and A = 200 for weak 
lensing analyses. 

The most direct probe of the cluster mass is weak grav¬ 
itational lensing (for a review, see Bartelmann & Schnei¬ 
der 2001; Refregier 2003; Schneider 2006b; Hoekstra & Jain 
2008; Massey et al. 2010), which is independent of the dy¬ 
namical state of the cluster and depends only on gravity. 
By statistically detecting correlations in the images of more 
distant galaxies whose light has passed the objects of inter¬ 
est, we learn about the gravitational fields of these objects, 
meaning weak lensing measurements are equally sensitive to 
dark matter and baryonic matter (Schneider 2006a). Galaxy 
clusters provided the first demonstration of weak lensing 
(Tyson et al. 1990) and measuring galaxy cluster masses, 
both on an individual basis and in stacks of similar objects, 
is an active field of research. However, weak lensing mea¬ 
surements have their own systematics, including projection 
effects and the effects of nonspherical haloes and cluster sub¬ 
structure (e.g., Corless & King 2007; Marian et al. 2010; 
Meneghetti et al. 2010; Becker & Kravtsov 2011; van Uitert 
et al. 2012). 

Several previous studies have used weak lensing to cal¬ 
ibrate easier to measure, but harder to interpret, cluster 
observables from multiple wavebands. We focus here on cal¬ 
ibrations of X-ray mass scalings, since these measurements 
have been used as inputs in recent broad cosmological mea¬ 
surements. See, e.g., the Arnaud et al. (2010) comparison 
used as input to the measurements of Planck (Planck Col¬ 
laboration et al. 2014b, a); the Arnaud et al. 2010 relation 
was also used to generate the X-ray scaled masses in the 
RBC catalogue. Unless otherwise noted, all quantities are 
reported at Ac = 500. We report bias here and throughout 


the paper as 

Mx/Mwl = 1 - 5 = g. (1) 

• Gruen et al. (2014) measure weak lensing masses for 
twelve SZ-selected galaxy clusters from the South Pole Tele¬ 
scope and compare them to results from the literature from 
Chandra, ROSAT, and XMM-Newton. The clusters range 
from redshifts of 0.1 to 0.7, with most between 0.3 and 0.4, 
and from masses of 0.4 — 2 x lO^^/ifg^M©. They find results 
for a fit bias consistent with no hydrostatic mass bias ( 21(^5 
per cent) but also marginally consistent with the expected 
~ 20% bias from simulation. 

• In an analysis of four relaxed clusters with joint Suzaku 

and Subaru observations, Okabe et al. (2014) compute hy¬ 
drostatic masses and find a ~ 24 ± 7% average bias at 
Ac = 500, increasing to ~ 60 ± 10% at the virial overden¬ 
sity. Redshifts range from 0.05 to 0.25 and masses (in this 
case M 2 oo,c) from 0.4 —1.7 x They, too, observe 

increased bias at larger radii. 

• Israel et al. (2014) examine 8 clusters from the 400d sur¬ 
vey. The clusters have redshifts 0.39 < 2 < 0.8 and masses 

1 — 6 X lO^^M©. They compare weak lensing masses with 
hydrostatic mass estimates and find results consistent with 
no bias, and with somewhat lower scatter than expected. 

• Donahue et al. (2014) investigate 25 CLASH clusters 
using weak and strong lensing profiles in conjunction with 
XMM-Newton and Chandra hydrostatic profiles combined 
as in Mahdavi et al. (2013). The CLASH clusters range from 
0.2 < 2 < 0.9 and masses M 25 qo,c from 1 — 10 x lO^^/i^i/M©. 
Donahue et al. find an average bias, but the magnitude and 
sign of the bias depends on the analysis type and the aper¬ 
ture used for the measurement. For the 0.5 Mpc measure¬ 
ment, the most similar to our measurements here, they find 
an average bias of 12 ± 7 per cent for weak lensing (consis¬ 
tently with previous measurements) but —11 ±7 for a strong 
and weak lensing combined measurement. The average bias 
increases to 20 to 30 per cent at 1 Mpc. 

• Applegate et al. (2015) compare twelve relaxed clus¬ 
ters with X-ray data from Chandra with weak lensing 
masses from the Weighing the Giants programme. The clus¬ 
ters have redshifts 0.15 — 0.54 and masses M 25 oo,c from 

2 — 12 X 10^'^M©. The average bias is found to be —4±8 per 
cent, with the weak lensing masses smaller than the X-ray 
masses but consistent within the errors. The same group 
made a measurement (von der Linden et al. 2014b) using 
this set of twelve clusters plus ten others and compared the 
weak lensing measurements to SZ masses instead of X-ray 
masses, finding a bias of 31 ± 7 per cent, with the X-ray 
masses larger than the weak lensing masses. 

Some earlier studies were performed with less constrain¬ 
ing power, such as Foex et al. (2012) from EXCPRES and 
Zhang et al. (2007, 2008, 2010) from LoCuSS. Comparisons 
are also performed against other X-ray measurements, such 
as temperature; for example, Israel et al. (2015) find a bias of 
about 35 per cent comparing to temperature, but note that 
this is different by 15-20 per cent than the bias measured 
with Chandra instead of XMM-Newton, pointing to the in¬ 
fluence of X-ray calibrations on the bias. Other studies (e.g. 
von der Linden et al. 2014b) have calibrated masses from 
cosmic microwave background data, which are ultimately 
scaled from X-ray masses, but are not a direct comparison. 
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Still others (e.g. Giles et al. 2015; Mahdavi et al. 2013) fit 
the X-ray scalings directly with weak lensing data, without 
first passing through a theoretically motivated proxy-mass 
relationship. 

In this paper, we use weak lensing to measure 
the masses of a subsample of the MCXC X-ray clus¬ 
ter catalogue (Piffaretti et al. 2011) called the RE- 
FLEX-hBCS-hCIZA (RBC) catalogue (e.g., Hajian et al. 
2013), which has a well-defined selection function with a 
set of predicted masses from scaled X-ray luminosities (Ar- 
naud et al. 2010). (We note that due to extinction cuts for 
lensing quality we do not include any GIZA clusters, but we 
retain the acronym since the catalogue has a common for¬ 
mat and standardized measurements.) We recalibrate these 
masses using the subset of RBC clusters that are found in 
the region of the Sloan Digital Sky Survey (SDSS; York et al. 
2000), using the weak lensing signal measured using back¬ 
ground galaxy shapes and photometric redshifts from SDSS. 
We also have a set of hydrodynamic simulations with similar 
characteristics to the RBC catalogue, which we will exploit 
to produce a model for the weak lensing signal including 
non-trivially complex physical effects such as triaxiality and 
substructure. As the RBC is an all-sky survey, and thus over¬ 
laps with many upcoming surveys and smaller projects, we 
hope this calibration will be useful for many groups. 

We cover the theoretical background of weak lensing 
in Section 2. Section 3 describes our data and simulations, 
and Section 4 shows the results of our measurements. We 
summarize the mass calibration in Section 5. 


2 METHODOLOGY 
2.1 Weak lensing 

Weak gravitational lensing leads to a correlation between 
the shapes of galaxies due to the gravitational influence of 
the matter between those objects and the observer. A con¬ 
centrated overdensity of matter, such as a cluster, will alter 
the apparent shape of the background galaxies to produce 
tangential correlations: due to the gravitational field of the 
overdensity, the average galaxy shape in an annulus will ap¬ 
pear to lie preferentially along the annulus itself, rather than 
randomly or along the radial direction. We call this distor¬ 
tion of the shape a shear. Weak gravitational lensing is called 
“weak” because these shears are small. 

We work in the thin-lens limit, where the line-of-sight 
extent of the lens is much smaller than both the distance 
between the observer and the lens and the distance between 
the lens and the background galaxy. A lens with surface 
mass density S(R) at redshift zi will induce on a background 
galaxy at redshift z a tangential shear 7 t (Schneider 2006a): 


WR) 


S(< R) - S(R) 

Ecr(zz, 2:) 


( 2 ) 


The shear is proportional to the difference between the aver¬ 
age surface mass density interior to radius R (E(< R)) and 
the average surface mass density at the radius R (E(R)). No 
assumption of spherical symmetry is necessary: the equation 
is exact for any overdensity (or underdensity) of matter as 
long as the average around the entire annulus or disk is con¬ 
sidered. The shear is also inversely proportional to the crit¬ 


ical surface mass density, or Ecr, which is given by a scaled 
combination of angular diameter distances: 


Ecr — 


Da (2 


47rG (1 - I - ziWDa{zi)Da[zi, z) 


(3) 


where the extra factor of (1 -|- Zi)~^ comes from our use 
of comoving coordinates. A more formal definition of weak 
lensing, as opposed to strong lensing, is the regime where 
E(R)/Ecr(R) < 1. 

We combine the noisy per-object galaxy shape measure¬ 
ments into a single statistic via an estimator 


AE = C{R) 




(4) 


The optimal weighting for this equation includes both the 
per-object shape noise weighting uiqshape and the critical 
surface mass density: 

~ rUi,shape (5) 


(Johnston et al. 2007; George et al. 2012). The critical sur¬ 
face mass density optimally reduces the scatter caused by 
the different lensing geometries of background galaxies at 
different redshifts. The per-object shape noise weighting will 
be described in Section 3.2. 

The optimal weighting also includes what is known as 
a boost factor correction, C{R). This accounts for the con¬ 
tamination of the background galaxy sample by galaxies as¬ 
sociated with the clusters themselves, which have a shear 
expectation value of 0, rather than the positive value ex¬ 
pected for lensed galaxies. Usually this correction is made 
via a comparison to weighted number densities (using the 
weights in Eq. 5) around random points distributed in the 
same way as the lenses (Sheldon et al. 2004). 

For our purposes, with such a small number of lenses 
covering the full extent of the somewhat heterogeneous back¬ 
ground galaxy sample, random points would likely not pro¬ 
duce the expected distribution of background galaxy densi¬ 
ties. Instead, in this analysis we take the average density in 
a large-scale annulus around the clusters - from 18 to 
Mpc - which is larger than our field of interest and unlikely 
to be contaminated by cluster galaxies, but is still physi¬ 
cally close enough on the sky to probe the same field number 
densities as the lenses themselves. Since survey edges inter¬ 
sect this largest annulus (and, to some degree, the annuli 
closer to the cluster centers), we cannot use the expected 
area of the annulus to modify the galaxy densities. Rather, 
we measure the number of red galaxies with photometric 
redshift above 0.5. For red galaxies at higher redshift, our 
photometric redshift estimates are secure and contamina¬ 
tion from cluster member galaxies is negligible. This sample 
is too small to detect the lensing effect above the noise, but 
it is large enough to be suitable for measuring the source 
number density. For a constant density of red galaxies, a 
fair assumption for the SDSS catalogue, the area is simply 
proportional to the number count of red galaxies in the an¬ 
nulus. This gives us a boost factor 


C{Ri <R<R2) 


f Yred(18 < Rj < 40) \ 

\Nred{Rl < R < R 2 ) / Ej;lg<flj<40 

( 6 ) 


Because of the extreme contamination due to survey 
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Figure 1. Sky coverage of the SDSS shape catalogue (in black) 
and the overlapping RBC clusters (red Xs). 

edges in the sparsely-covered southern Galactic cap region 
of SDSS (the three stripes broken at a right ascension of 0 in 
Fig. 1), even this method was insufficient to accurateiy mea¬ 
sure the boost factor in that region. The boost factor should 
be consistent for clusters of the same richness and redshift, 
so we should be able to use the 150 clusters of the northern 
cap to stand in for the 16 southern cap clusters if the rich¬ 
ness and redshift distributions are consistent. We find that 
the richness distribution is consistent, but the redshift dis¬ 
tribution shows a larger concentration of low-redshift clus¬ 
ters compared to the northern cap. We at first attempted 
measuring the boost factor in redshift slices and reweight¬ 
ing to obtain the global boost factor, but we found that 
the reweighted boost factors had approximately the same 
errors relative to the all-redshift measurement as a simple 
unreweighted combination of the redshift slices, so we use 
the directly measured boost factor for clusters of all red- 
shifts in the northern cap as the boost factor for the entire 
sample. This result holds when we bin the clusters in X-ray 
mass. 

Errors on the signal are estimated as the standard devi¬ 
ation of the mean for each sum in Equations 4 and 6. Since 
we are looking here at small scales around widely-separated 
sources, we do not expect a noticeable covariance between 
adjacent radial bins. 

The AE measurements as a function of radius give us 
an estimate of the 2-dimensional surface mass distribution 
in the region of the cluster, which we must turn into a 
3-dimensional mass contained within some region of inter¬ 
est. There are nonparametric methods of performing this 
inversion, but since we believe we know the density pro¬ 
files of galaxy clusters, interpretation is made simpler by 
using preexisting 3-dimensional models projected down to 
2-dimensional surface mass densities, which we can then 
compare to our measured AS profiles. In this work, we will 
compare our results with theoretical NFW profiles and with 


results from hydrodynamic simulations described in Section 
3.3. 


2.2 Miscentering 

We measure the weak lensing signal relative to the reported 
galaxy cluster center, which is the X-ray centroid from the 
X-ray catalogue described in Section 3.1. Errors in the deter¬ 
mination of the cluster center will suppress the weak lensing 
signal on small scales. Since we are taking averages of shears 
in annuli, an annulus around a location other than the peak 
will contain a small region of higher shear than expected and 
a much larger region of lower shear than expected. This ef¬ 
fect reduces the average shear in the annulus, and the effect 
size is greater at smaller radii, where the shear is chang¬ 
ing rapidly with radius. If not accounted for directly, this 
suppression will cause errors in our final mass measurement 
(Johnston et al. 2007; Mandelbaum et al. 2008a; Rozo et al. 
2011) that can be tens of per cent, comparable to the size 
of the mass bias that we are trying to measure. 

Errors in centroid determination can come from phys¬ 
ical sources (unrelaxed clusters whose center of mass may 
be physically separated from the X-ray peak) and from 
measurement uncertainty (the finite resolution of the X-ray 
map). For optical cluster catalogues, some fraction of the 
clusters are exactly centered, due to the exact centering of 
some brightest cluster galaxies. The rest have a distribution 
of Rs, the distance between the true center and the observed 
center, usually written as a Rayleigh distribution (Johnston 
et al. 2007) 

P(R.) = ^exp(^-i(R./a,)^^ (7) 

which is the average absolute offset for a measurement with 
Gaussian error in both the x and y directions. For the X- 
ray measurements, we have no perfectly centered clusters 
since centroid measurement uncertainty affects them all, so 
we may use P{Rs) directly for our sample as a whole. For 
as, the width of the distribution of centroiding errors, we 
use both 0.25' and 0.5', related to the average resolution of 
the X-ray maps used in the creation of the RBC catalogue 
described in section 3.1; we also tested several other val¬ 
ues, as discussed in section 4, but these values encompass 
the best fits to the data, as described further in section 4. 
Our theoretical mass profiles E(P) are effectively convolved 
with this distribution of centroid offsets before calculating 
the theoretical lensing signals, AE. For our simulation re¬ 
sults, this convolution is done on a cluster-by-cluster basis 
using the simulation redshifts, since we must convert from 
angular to physical coordinates. For the NFW profiles, due 
to computation time limits, the clusters are grouped into 
redshift slices containing 10-20 clusters, sometimes less at 
the high-redshift end of the sample. The average weight and 
angular diameter distance of the groups are calculated us¬ 
ing weights based on the lensing geometry, and a weighted 
summed profile is constructed using the results from each 
group of clusters instead of each cluster individually. 
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3 DATA AND SIMULATIONS 
3.1 The RBC cluster sample 

The cluster catalogue that we used is the RBC catalogue, 
which is a sub-sample of the MCXC X-ray cluster cat¬ 
alogue (Piffaretti et al. 2011). This sub-sample combines 
the REFLEX (Bohringer et al. 2007), BCS (Ebeling et al. 
1998, 2000) and CIZA (Ebeling et al. 2002; Kocevski et al. 
2007) flux limited catalogues from the RoSAT All Sky Sur¬ 
vey (RASS; Voges et al. 1999). (We do not include any 
CIZA clusters after our lensing cuts-see below for more 
information.) We assume the RBC catalogue is complete 
above a flux limit of ~ 3 x 10“^^/ergs“^cm“^, or 4.4 x 
10“^^/erg s“^cm~^ for the BC subsample, so we look at clus¬ 
ters only above this flux limit. These are the same cuts used 
in Mantz et al. (2010), which also contains more informa¬ 
tion on possible incompletenesses. We consider only the 166 
clusters that overlap with the SDSS lensing footprint. The 
cluster masses are taken from the MCXC catalogue, which 
uses a Lx — M scaling relation (Malmquist bias corrected) 
from the REXCESS cluster sample (Pratt et al. 2009). The 
scaling relation is 


h{z) 


- 7/3 /_ 


\10"‘'^ergs“^ 


= C 


Ms, 


2 X lOi^Mo 


( 8 ) 


where h{z) = [nM(l -I- z)^ -f and log(C') = 0.274. 

Additionally, there is scatter about this relation of ~ 24% 
in mass (or ~ 38 % in Lx Pratt et al. 2009). The scaled 
cluster masses in the Pratt et al. (2009) relation were ini¬ 
tially calculated under the assumption of hydrostatic equi¬ 
librium (HSE Pointecouteau et al. 2005) for a sample of 
10 relaxed clusters. They were then used to calibrate the 
Tx — M (Kravtsov et al. 2006) relation in Arnaud et al. 
(2007) and updated in Arnaud et al. (2010). This larger set 
of clusters with a lx — M relation were then used to infer 
a relationship between M and Lx. This scaling relation was 
then extended to an even larger sample of clusters in the 
analysis of Pratt et al. (2009). Therefore, the assumption of 
HSE is still present in the masses presented in the MCXC 
catalogue, albeit at some remove. We do not account for an 
HSE bias or any additional bias that would result from using 
this Lx — M relation, since that bias is what we are trying 
to measure. 

The cluster sample in the SDSS region ranges from 
2 ~ 0.04 — 0.4 with a mean redshift of 0.14. The re¬ 
ported masses from scaled luminosities range from M 500 ~ 
0.6 — 10.5 X Mq, with a mean of 3.4 x Mq. 

A map of the SDSS coverage area and the overlapping RBC 
clusters is shown in Fig. 1, and the catalogue’s redshift and 
luminosity distribution are shown in Fig. 2. We note that no 
CIZA clusters appear in the masked catalogue, since they 
are located in regions of high extinction where photometric 
redshifts cannot be reliably determined. 


3.2 Weak lensing data 

3.2.1 Galaxy shape sample 

We use a shape catalogue measured using the re- 
Gaussianization software (Hirata & Seljak 2003) applied to 
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Figure 2. Luminosity and redshift coverage of the RBC clusters 
that overlap with the SDSS source catalogue used in this paper. 
The horizontal line indicates the approximate location of the split 
between our two mass bins. 


data from the Sloan Digital Sky Survey (SDSS) data re¬ 
lease 8 (York et al. 2000; Aihara et al. 2011). This cata¬ 
logue is described in detail in Reyes et al. (2012), includ¬ 
ing details of shear calibration and error estimation. The 
software calculates adaptive second-order moments for the 
galaxy images and PSFs by effectively fitting them to ellip¬ 
tical Gaussians, and those moments are then combined to 
generate PSF-corrected distortion measurements including 
a correction for low-order non-Gaussianity of the galaxies 
and PSFs. We turn the distortions into shears by dividing 
by a shear responsivity factor TZ (Bernstein & Jarvis 2002) 
which for our shear catalogue is approximately 1.8 (Reyes 
et al. 2012). This factor describes the average sensitivity 
of the shape sample to an applied shear, as shears do not 
add linearly (Bernstein & Jarvis 2002). After cuts removing 
badly measured or unmeasurable objects, the catalogue con¬ 
tains 39 million galaxies in an area of approximately 9,000 
square degrees, for an overall unweighted number density 
of 1.2 galaxies/arcmin^. This shear catalogue was tested ex¬ 
tensively in Reyes et al. (2012). Additionally, further testing 
was done using simulated images, which found that com¬ 
bined errors from PSF correction, noise bias, and selection 
bias gave a systematic shape error budget of 3.5 per cent 
(Mandelbaum et al. 2012); combined with other calibration 
errors such as photo-z error and stellar contamination, AE 
profiles calculated using this shape sample should have a to¬ 
tal systematic error budget of 4 per cent including all shear 
and photo -2 biases (Mandelbaum et al. 2013), below the 
statistical error of our measurement. 

The galaxies have a per-object shape weight of 

Wi,shape = , 2 ' (9) 

^rms “r ^ e,i 

which is a combination of the individual measurement error 
(Je,i and the intrinsic dispersion of galaxy shapes Crms, which 
is 0.36 for this sample. 
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3.2.2 Galaxy redshifts 

The galaxy photometric redshifts are presented in Naka- 
jima et al. (2012). The measurements were made with the 
Zurich Extragalactic Bayesian Redshift Analyzer, or ZE¬ 
BRA, a template-fitting software (Feldmann et al. 2006). 
Thirty-one templates were used, with six templates (four 
observed SEDs and two synthetic blue galaxy spectra) from 
Benitez (2000) and twenty-five additional templates created 
via interpolation between neighboring pairs of the six start¬ 
ing templates. Starburst-type galaxies, representing 10 per 
cent of the sample, were removed due to unreliability after 
comparison against spectroscopic and high-quality photo¬ 
metric redshifts. The remaining galaxies have a known bias 
for 2 > 0.4, approximately the mean redshift of the cat¬ 
alogue. The bias in the lensing signal due to background 
galaxy photo -2 bias and scatter can be corrected since we 
know the true dN/dz (Nakajima et al. 2012). For this cluster 
sample, in which the majority of galaxy cluster lenses are at 
2 < 0.1, the typical correction is small, from 0 to 5 per cent 
depending on binning scheme. 

The noise in the photo -2 estimator also makes our opti¬ 
mal weighting, Eq. (5), somewhat less optimal, correspond¬ 
ing to an approximately 10 per cent increase in the error 
bars compared to a case with no noise in the source galaxy 
redshifts (Nakajima et al. 2012; Mandelbaum et al. 2008b). 
However, this is still a factor of ~ 2 better than the case 
without optimal weighting. 

3.3 Cosmological Hydrodynamic Simulations 

We used simulations of galaxy clusters to interpret the weak 
lensing observations. The clusters were extracted from sim¬ 
ulated cosmological volumes (L = 165 Mpc/h) such that 
there was a sufficient statistical sample of clusters to cal¬ 
culate the lensing signal. These simulations were run with 
a modified version of the GADGET-2 smoothed particle 
hydrodynamics (SPH) code (Springel 2005). This version 
of the GADGET-2 code included sub-grid models for ac¬ 
tive galactic nuclei (AGN) feedback (for more details see 
Battaglia et al. 2010), radiative cooling, star formation, 
galactic winds, supernova feedback (for more details see 
Springel & Hernquist 2003), and cosmic ray physics (for 
more details see Pfrommer et al. 2006; Enfilin et al. 2007; 
Jubelgas et al. 2008). We ran 10 165 Mpc//i boxes with 
a resolution of 256^ gas and dark matter (DM) particles 
which yields a mass resolution of Mgas = 3.2 x IO^Mq/Zi 
and Mdm = 1.54 x IO^^Mq/Zi. The cosmological parame¬ 
ters used for these simulations were Hm = floM + Hb = 0.25, 
Hb = 0.043, Ha = 0.75, i7o = 100/ikms"^ Mpc“\ h = 0.72, 
Us = 0.96 and ng = 0.8. 

At each redshift snapshot the halos were identified and 
their properties were calculated in two steps: 

• The halos were found using a friends of friends algo¬ 
rithm (Huchra & Geller 1982). 

• For each halo the center of mass was computed itera¬ 
tively, and then the spherical overdensity mass (Ma) and 
radius (Ra), where A refers to the multiplicative factor ap¬ 
plied to the critical density, pcr(z) = 3JT|[Hm(1 + z)^ + 
HA]/(87rG). 

For each redshift slice we project the mass distributions of 


the halos. These projections are 8x8 Mpc/h comoving cen¬ 
tered on the cluster’s COM. We made separate gas, stellar, 
and DM components projected down the entire 165 Mpc//i 
length. Then we summed the projections of all three com¬ 
ponents for our analyses. 

The simulated stacked AE profiles that we compare to 
the RBC weak lensing observations are calculated in the 
following way. We apply a selection function to the simu¬ 
lated halos which uses the scaled masses and redshifts from 
the RBC catalogue, since the clusters in the RBC sample 
are flux limited. We model this selection using the initial 
spherical overdensity masses from the simulations and then 
we apply a Monte Carlo method that scales each of these 
masses by a randomly selected number drawn from a Gaus¬ 
sian distribution with a mean of 1 and a dispersion that 
matches the Lx — M scatter in mass of 24%. This correc¬ 
tion models the Eddington bias resulting from the scatter 
in the scaling relation used to calculate the masses in the 
RBC catalogue. These modeled observed masses are used to 
determine which halos are placed in a given mass bin. We 
then apply the redshift cut, which is well determined for a 
given mass from the RBC selection function. We emphasize 
that we do not ever directly predict luminosity quantities: 
we always work in terms of mass for these predictions. 

For each halo we compute the AE profile with miscen- 
tering included. The probability distribution for miscenter- 
ing is described by Eq. 7. In the simulations this miscenter- 
ing is modeled by a Monte Carlo method. For each cluster 
we displaced the center of the cluster by 

Rs = (10) 

where a;rand is a random draw from a uniform distribution. 
Eq. 10 results from integrating Eq. 7 to get the cumulative 
distribution function and solving for Ra. We used values of 
0.25' and 0.5' for as, and the motivation for this choice is 
discussed in Sec. 4. For each halo we found a new center 
10 times and then calculated the average AE profile given 
these 10 new centers. 

The final simulated AE profiles for a given mass were 
calculated by averaging over all halos within a given mass 
range. This includes halos at different redshift outputs in the 
simulations. We calculate a weighted final average across the 
redshift outputs, which accounts for both the weak lensing 
geometry-related weights and the volume factor (comoving 
distance squared) associated with including multiple red¬ 
shifts. 

3.3.1 Eddington and Malmquist bias 

When measuring the Lx — M relation in real data, the initial 
catalogue corrected for the Malmquist bias of their sample as 
follows. The calibration was measured on an initial sample, 
and then the relation was altered using simulation to repre¬ 
sent what the calibration would have been if no Malmquist 
bias was present in the observation itself. This procedure 
yields a relation which gives a correct mean mass estimate 
for each Lx measurement. 

However, once we start to bin in the masses generated 
from the Lx — M relation, we introduce Eddington bias due 
to the scatter in each individual mass. We model this ef¬ 
fect by adding a random number drawn from a distribution 
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with the expected scatter to the known true masses of our 
simulated clusters before binning, as described above; this 
mimics the effect of binning the scattered (but no longer 
Malmquist-biased) masses of the original catalogue. 


4 RESULTS AND DISCUSSION 

We show results for the cluster sample as a whole and also 
for two mass bins, divided as shown in Table 1. The cluster 
sample as a whole has the smallest statistical error, since it 
contains the most objects, with a total of 1.2 million lens- 
source pairs; however, since it spans more than an order 
of magnitude in mass, interpretation of results such as a 
mean mass can be difficult. But making our bins too narrow 
will also result in systematic errors, since we are relying 
on the advantage of stacking many objects to reduce the 
effects of stochastic uncorrelated large-scale structure and 
triaxiality. Even in the two-bin case, our higher-mass bin 
has fewer than the 100 clusters suggested by Corless & King 
(2009) to effectively reduce triaxial bias, and a similar effect 
is likely to come into play for other complications whose 
impact we expect to be reduced by stacking, such as large- 
scale structure. Therefore, we do not believe results from 
smaller bins will be reliable. As stated previously, the sample 
as a whole is a flux-selected sample (in the X-ray flux), which 
means that the mean mass is higher at high redshift than 
at low redshift. However, when we split into two mass bins, 
the higher mass bin has a high enough mass cut that it is 
more similar to a volume-limited sample than a flux-limited 
sample, as shown in Fig. 2. We also note that the lower mass 
bin has no clusters at the higher end of our redshift range. 

In Fig. 3, we show the weak lensing signal around the 
RBC clusters in one or two bins split by X-ray mass as de¬ 
scribed in Table 1, in addition to the expected signal from 
the simulations if there is no mass bias {q = 1) and if there is 
a 25 per cent mass bias {q = 0.75). We expect a slight decre¬ 
ment between the data and the model due to the mass bias 
(from assumptions of HSE as well as any calibration prob¬ 
lems or other sources of error), assuming the cosmology in 
the simulations is correct; the predictions are low compared 
to the data, though the difference is not large enough to be 
conclusive by eye. We show both the 0.25' miscentering off¬ 
set and 0.5' offset, since both provide reasonable agreement 
with the data, as we discuss later in this section. The ratios 
between the various predictions or data and the 0.25' off¬ 
set case are also shown for Bin 1/1 to facilitate comparison. 
The value for the q = 0.75 case in Bin 1/1 is 22.2, com¬ 
pared with 28.9 for g = 1 with the same miscentering offset, 
showing the clear need for such a calibration assuming the 
simulations have the correct cosmology. There is no fitting 
involved in this plot: this is a simple prediction from clus¬ 
ters selected from the simulation using the same selection 
function as the data. 

Next, we fit the weak lensing signal with a simple 
spherically-symmetric NEW profile, miscentered using the 
distribution in Eq. 7 in redshift slices as described in Sec¬ 
tion 2.2. We use the radius range 0.1 — 3 Mpc for the 
NEW haloes; this outer range corresponds to approximately 
the virial radius for a cluster at the mean mass and red¬ 
shift, so we should have minimal contribution from large- 
scale structure. We show results from two different options 


for the width of the miscentering distribution, 0.25' and 0.5'. 
Both fit the data relatively well, with some bins and fits pre¬ 
ferring the lower offset and some the higher. Since we are 
eliding other possible sources of error as well (such as possi¬ 
ble insufficient modeling of the brightest cluster galaxy), we 
believe that these results represent the allowable range given 
our knowledge of these systems. Offsets of 1' were tried and 
ruled out by increased values, and since 0.25' was already 
too narrow for some bins we did not try lower widths. Then 
we fix the mass-concentration relation to the z = 0 relation 
from Bhattacharya et al. (2013), which gives us a typical 
concentration of C 200 ~ 4 — 4.5 for this mass range. The 
best-fitting NEW profile, and the range of fits comprising 
the ±1(T masses, are shown in Fig. 4. Parameters of the fit 
are shown in Table 1. 

Finally, we fit mass profiles based on our simulations, 
using a radius range of 0.1 — 5.5 Mpc, which is about 
twice the virial radius. Since our simulations include large- 
scale structure, we use the increased radial range for in¬ 
creased statistical power; the actual mass values we find dif¬ 
fer by only a fraction of the Icr error bars compared to the 
values we get for an upper radius limit of 3 Mpc, so we 
are confident this does not bias our final results. We produce 
average AE profiles for a set of mass bins extending 0.60 in 
base-ten log mass space, with endpoints spaced every 0.24 
in base-ten log space, to constrnct a set of overlapping mass 
bins containing a distribution of cluster masses. We then lin¬ 
early interpolate between each neighboring, overlapping pair 
of mass bins to produce a AE profile at intermediate masses, 
and spline interpolate in the radial direction to produce a 
smooth function. We fit for the interpolation parameters: 
the index of the lower bin n and the fraction / needed to 
produce the best-fit AE, such that 

AEdata — (1 — f)^^n + /AEn-t-1 (H) 

We then interpret the mean mass of the fit to be the cor¬ 
responding linear combination of the mean masses of the 
clusters in mass bins n and n -|- 1. 

In principle, as the AE profiles are a power of the 
mass, our linear combination should properly be a logarith¬ 
mic mean. However, in practice, our bins are finely spaced 
enough, and the curves noisy enough, that we can interpo¬ 
late linearly without affecting the precision of our results. 
We tested this by generating mass bins from the simula¬ 
tions with the same size as the bins used for fitting the 
mass, but with offsets from the endpoints of those mass 
bins. With both linear and logarithmic interpolation, we 
find that at low to intermediate masses our fitting rou¬ 
tine returns the expected answer. The relation breaks down 
above ~ 7 x Mq, where our fitting code begins to 

return systematically low masses. The discrepancy rapidly 
increases with true mass, reaching a bias of 20 per cent 
or greater at masses greater than ~ 8 x 10 ^^/i/q^Mq. This 
breakdown above ~ 7 x 10^'^Mq corresponds to approx¬ 
imately the most massive 10 per cent of our sample (10- 
15 objects, depending on where we draw the line between 
stochastic noise and systematic bias). This is close to the 
values found for the higher-mass bin of the sample when it 
is split into two bins, and is possibly a source of error for 
those points, but we do not expect this to be a dominant 
source of error even for this bin (an expectation held up by 
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Mass 

Average 

MB-corrected 

Centroid 

Simulation-fit 


Mass 

NFW 


Mass 

Bin 

Afcl 

range 

Mx 

Mx 

offset 

Mwl 

X^/dof 

ratio 

Mwl 

X^/dof 

ratio 

1/1 

166 

0.6-15 

2.8 ±0.1 

2.5 ±0.1 

0.5' 

o q-|-0.5 
'^•^-0.4 

20.4/19 

0 79+® 09 

o- ' ^-0.12 

3.0 ±0.3 

23.2/16 

0.83 ±0.09 






0.25' 

o ^^-|-0.3 
'^•^-0.4 

17.7/19 


2.9 ±0.4 

26.9/16 

0.88 ±0.11 

1/2 

116 

0.6-4 

2.1 ± 0.1 

1.9± 0.1 

0.5' 

2.3 ±0.2 

16.3/19 

n S4+0'll 

U »4_q Qg 

2.3 ±0.2 

13.1/16 

0.84 ±0.09 






0.25' 

2.4 ±0.2 

23.9/19 

0.82tJ;“7 

2.2 ± 0.2 

14.1/16 

0.88 ±0.10 

2/2 

50 

4-15 

5.6 ±0.2 

4.8 ±0.2 

0.5' 

6 7+1-0 
^•'-0.7 

21.3/19 

0.72 ± 0.08 

5.6 ± 1.0 

24.8/16 

0.87 ± 0.16 






0.25' 

c 4+0.9 

21.4/19 

O- ' 0_Q Q7 

5.0 ± 1.0 

30.6/16 

0.96 ±0.21 


Table 1. Binning scheme for the RBC clusters, with characteristics of the X-ray masses from scaled luminosities (both the raw average 
and the average once we have corrected for Malmquist bias), masses fit from the weak lensing data using two lensing models (theoretical 
NFW and results from hydrodynamic simulations) for two different miscentering parameters, goodness-of-fit statistics, and the resulting 
X-ray to weak lensing mass ratios. All masses are in units of Mq and are lensing-weighted mean masses, rather than true mean 

masses; our results for true mean masses are consistent, but require larger corrections to obtain, so we show only the more reliable 
lensing-weighted averages. 
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Figure 3. Weak lensing measurements (in red) and predicted signal from simulations (solid line for the 0.25^offset, dashed line for the 
O.S^offset, and dot-dashed for the 0.25^offset with a calibration of g = 0.75) for the RBC clusters in two binning schemes with bins split 
by X-ray mass from scaled luminosities. Errors are statistical errors on the mean signal for both data sets, using y/n scaling. The bottom 
right panel shows the ratio between the various curves and the 0.25'offset curve for Bin 1/1 only, to better highlight the changes caused 
by these effects. 
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Figure 4. Weak lensing measurements (in red), best-fit NFW models with Icr error bars (in green) and best-fit interpolated simulation 
fit with 1(T error bars (in blue) for the RBC clusters in one or two X-ray scaled mass bins. Weak lensing errors are statistical errors; the 
NFW and simulation-fit error regions are defined by the Icr mass errors from the fitting procedure. Parameters for the fits are shown in 
Table 1. 


the general consistency of the results with the other bins 
and fitting methods). 

The interpolated curves and their error region are 
shown in green in Fig. 5. The interpolated simulation fits 
have a lower (or approximately equal) yf /Aoi for almost 
all fits, except the less massive of the two bins in the two- 
bin case for one value of the miscentering offset. By eye, 
the NFW profile is clearly “peakier” than the interpolated 
simulation fits, which we expect since baryonic effects in the 
hydrodynamic simulation make the cluster mass profiles less 
cuspy (e.g. Duffy et al. 2010). The data appears more consis¬ 
tent with the simulation trend at small radii, though not to 
a statistically significant level. This effect is not strongly de¬ 
generate with the miscentering effect, which tends to move 
the peak in radius, but not broaden it sufficiently to over¬ 
come the discrepancy in the small-scale behavior of the den¬ 
sity distribution. We did not investigate altering the mass- 
concentration relation of our NFWs to better fit the simula¬ 
tion results, as the simulations look non-NFW-like in other 
ways as well. 

A number of corrections are necessary to go from the 


raw X-ray averages and the raw lensing data to the mass 
comparison that is our final result. First, the lensing sig¬ 
nal must be corrected for the known errors in our photo -2 
catalogue, as defined in section 3.2.2. After we fit our para¬ 
metric models (either simulation or analytic NFW) to the 
data, we obtain a “mass” which is somewhere between the 
median and the mean of the actual distribution of masses 
that went into the signal, since the lensing signals are not 
linearly proportional to the mass; we must define an addi¬ 
tional correction to take our fitted masses and turn them 
into real mean masses. For the NFW case, we simulated a 
signal made of perfect NFW haloes with the mass distribu¬ 
tion of the RBC clusters, with several different multiplica¬ 
tive biases to mimic the bias effect, and compared our fit 
values to the known mean mass to calibrate this effect sepa¬ 
rately from errors caused by profile differences. The results 
were very stable across bias values covering the range found 
in this work. The largest correction is for the single mass 
bin case, as expected, with an 11 per cent calibration fac¬ 
tor, while the two-bin case has 6 per cent and 2 per cent 
for the lower-mass and higher-mass bins respectively. For 
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Figure 5. Mass ratios from the cluster bins of Table 1 for NFW 
profile fits (green) and interpolated simulation predictions (blue), 
plotted against a weighted average of the X-ray luminosity mass 
for each bin, slightly offset in the x-direction for clarity. Small 
points are the results from the 2-bin case, large points from the 
1-bin case. Error bars are Itr statistical errors from the fitting 
procedure on the y-axis and scatter within the X-ray bin on the 
x-axis. The black line is unity. 


the simulation case, we fit the predictions shown in Fig. 3 
with our interpolation scheme and use the resulting bias 
as our correction, which is 6 per cent for bin 1/1, 0.1 per 
cent for bin 1/2, and 8.5 per cent for bin 2/2. Finally, we 
must weight the X-ray masses from the RBC catalogue by 
the expected geometric leasing factors (nearby clusters take 
up more area on the sky, meaning more lens-source pairs 
in the sums, and also higher lensing efficiency due to the 
distance ratios involved) and correct the resulting average 
for Malmquist bias. We use the relationship between true 
mass and “observed” (scattered) mass in our simulations to 
measure the Malmquist bias correction factor, while the ge¬ 
ometric lensing weights are purely analytic once we fix the 
cosmology to our hducial model. The corrections range from 
0 to 5 per cent for photo-z; from 0 to 11 per cent for lensing 
mass vs mean mass differences; and from 9 to 15 per cent 
for Malmquist bias. 

We note that, while both the X-ray and weak lensing 
masses are reported at Ac = 500, the weak lensing masses 
are the product of a parametric fit that goes to larger radius, 
so the masses are not measured on the same aperture. 

Finally, in Figure 5, we show the two mass determina¬ 
tions (NFW and simulation) for the two miscentering pa¬ 
rameters in comparison to the expected average mass de¬ 
rived from the X-ray luminosity masses. The simulation hts 
and NFW fits show a clear trend for a larger weak lensing 
mass than X-ray mass, with a size of approximately 20 per 
cent, consistent with our expectations from simulations. 


5 CONCLUSION 

We have used weak lensing measurements to calibrate the 
average masses for 166 clusters with X-ray measurements in 
the RBC catalogue. We find evidence that the X-ray masses 
from scaled luminosities in the RBC catalogue are approxi¬ 
mately 15 to 30 per cent lower than the weak lensing masses 
over the range of masses probed in this cluster sample (0.6 to 
10.5 X Mq), with slight but not statistically signif¬ 

icant evidence of a mass-dependent bias for the simulation 
fits, which we consider more reliable. We cannot say with 
certainty what fraction of this mass bias is due to the as¬ 
sumption of hydrostatic equilibrium: given the discrepancy 
in X-ray masses between different analysis methods (Rozo 
et al. 2014), it is likely that some of the bias is due to effects 
other than the assumption of hydrostatic equilibrium, such 
as differing X-ray calibrations, errors in the original Lx — M 
fit in addition to problems of hydrostatic equilibrium, and 
discrepancies between X-ray and weak lensing estimates of 
Rsoo- However, our analysis is agnostic to the sources of dif¬ 
ference between the X-ray and weak lensing estimates, and 
merely produces a corrected version of the masses given by 
scaled luminosity in he RBC catalogue. This measurement 
is in agreement with the expectations from simulations (e.g. 
Okabe et al. 2014; Donahue et al. 2014; Mantz et al. 2014). 

To estimate the basic implications for a cosmological 
measurement using cluster abundances, we use the mass 
function from Tinker et al. (2008) at z = 0 and the cosmo¬ 
logical parameters from Planck Collaboration et al. (2014a). 
Assuming use of all clusters down to a (real) mass limit of 
Mq, then if a mass bias corresponding to g = 0.78 
is ignored, a single-parameter cosmological constraint on erg 
would be biased low by 12 per cent. For higher masses, this 
rapidly becomes worse due to the steepening of the mass 
function. Obviously a real analysis would not vary only a 
single parameter, but this gives a rough indication of the 
magnitude of the problem. Since many cosmological analy¬ 
ses have smaller statistical errors than this, it will be im¬ 
portant to apply corrections for this mass bias rather than 
using X-ray-inferred masses directly. We defer a full analysis 
of the effect on multi-parameter cosmological constraints to 
future work. 

There are cosmological implications for this measure¬ 
ment with respect the previous Planck cluster cosmology 
analysis. In Planck Collaboration et al. (2013) the assumed 
value for the mass bias value is consistent with our result. 
The masses in Planck Collaboration et al. (2013) are cali¬ 
brated by the REXCESS X-ray masses (Planck Collabora¬ 
tion et al. 2011) just like the MCXC and RBC catalogues. 
Our results are inconsistent with the 40-50 per cent mass 
bias required to bring the cosmological parameters from 
the Planck cluster analysis into agreement with the primary 
CMB, and agree well with previous analyses such as von der 
Linden et al. (2014b) (30 per cent bias) and Hoekstra et al. 
(2015) (25 per cent bias). Our analysis does not rule out 
contributions for other systematics in the Y-M relationship 
used in Planck Collaboration et al. (2013). 

Our reported hydrostatic mass bias of 22/)®(, per cent is 
within 2(7 of the recent results from the LoCuSS subsample 
presented in Smith et al. (2015) from mass measurements in 
Okabe & Smith (2015). They find a 5 ± 5 per cent hydro¬ 
static mass bias. When comparing these results one should 
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keep two things in mind: the X-ray masses use different scal¬ 
ing relations and our results and the results of Okabe & 
Smith 2015 apply to different cluster samples. In relation 
to the Planck cluster cosmology results, the bias we calcu¬ 
late needs to be carefully propagated through the previous 
Planck calibrations of cluster masses and the Smith et al. 
(2015) bias should be carefully corrected for the Eddington 
bias inherent in the Planck cluster masses (Battaglia et al. 
2015). 

Beyond the Planck cluster analysis our results are an 
important calibration for any results that use the RBC clus¬ 
ter mass estimates in their cosmological parameter estima¬ 
tion. We have quantified a systematic bias that other pre¬ 
vious analyses (e.g. Hajian et al. 2013) had to marginalize 
over with broad priors. New or reanalyses of measurements 
that use the RBC cluster catalogue will signihcantly improve 
with our mass calibration. 
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